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ABSTRACT 

We compare the rms contrast in observed speckle reconstructed G-band images 
with synthetic filtergrams computed from two magneto-hydrodynamic simulation 
snapshots. The observations consist of 103 bursts of 80 frames each taken at the Dunn 
Solar Telescope (DST), sampled at twice the diffraction limit of the telescope. The 
speckle reconstructions account for the performance of the Adaptive Optics (AO) sys- 
tem at the DST to supply reliable photometry. We find a considerable discrepancy in 
the observed rms contrast of 14.1% for the best reconstructed images, and the syn- 
thetic rms contrast of 21.5% in a simulation snapshot thought to be representative of 
the quiet Sun. The areas of features in the synthetic filtergrams that have positive or 
negative contrast beyond the minimum and maximum values in the reconstructed im- 
ages have spatial scales that should be resolved. This leads us to conclude that there 
are fundamental differences in the rms G-band contrast between observed and com- 
puted filtergrams. On the basis of the substantially reduced granular contrast of 16.3% 
in the synthetic plage filtergram we speculate that the quiet-Sun may contain more 
weak magnetic field than previously thought. 

Subject headings: Sun: granulation — Sun: photosphere — techniques: high angular 
resolution — (magnetohydrodynamics:) MHD 



1. Introduction 

Magneto-hydrodynamic simulations are increasingly successful in reproducing the observable 
effects of magneto-convection on the solar surface. In achieving this goal they provide not only 
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insight into the physics that drives convective patterns at the surface, but also into the dynamics 
of the solar interior, where at best only helioseismological observations can probe. Although the 
simulations are typically claimed to be parameter free, and therefore cannot be adjusted at will to 
match observations, they include arbitrary elements in the choice of the physical approximations 
that are implemented to keep the numerical problem tractable. Significant limitations are restrict- 
ing radiative transfer to a small number of representative frequency bins and discrete directions, 
and strongly reducing th e Reynolds numbers in the simulation to values that are orders of magni - 
tude less than in reality Jstein & Nordlundlll998l : Ivogler et allboosi Ischaffenberger et allbood) . 
The latter restriction limits the smallest spatial scales at which structures can occur and how the 
dynamics on these scales can possibly feed back to larger scales. 

One of the major successes of solar convection simulations is the very accurate reproduction 
of spatially and tempor ally averaged shapes of weak photospheric lines (lAsplund et al.l l2000bl : 
Stein & Nordlundll2000l) . These line shapes, characterized by an asymmetry towards the blue, are 
the intricate product of correlation between up- and downflows, and their respective brightness 
differences and area asymmetries. Because the comparison is made on the basis of spatially and 
temporally averaged profiles spatial resolution in the observations is irrelevant. This is not true for 
the numerical resoluti on, which needs to be sufficiently refined to provide a good match to observed 
average line profiles (|Asplund et al.ll2000ah . However, the averaging may mask discrepancies in 
separate quantities, like brightness contrast, velocity variations, and size spectrum of the convective 
structures that could compensate each other in the averaging process. 

In concert with developments in theoretical modeling, a concomitant progress in observational 
techniques of the solar atmosphere has taken place. In particular, the advent of high order, high 
temporal bandwidth adaptive optics and post facto image reconstruction techniques, like speckle 
interferometry (e.g.Jvon der Luehelll993b and multi-object multi-frame blind deconvolution (e.g., 
van Noort et al.ll2005b . have greatly advanced our capabilities of imaging the solar surface with 
high fidelity. Together these developments should allow for a promising and more realistic com- 
parison between theoretical predictions and observations. Here we attempt to take advantage of 
developments in instrumentation in order to compare simulations and observations on the basis 
of only the rms contrast, in particular that in the G band. This eliminates the possible confusion 
of mixing velocity signatures with intensity contrast that is inherent in the comparison of spa- 
tially and temporally averaged line profile shapes, but of course has the drawback that we have to 
compare with spatially resolved observations. Thus we have to compensate for the image degra- 
dation imposed by the Earth's atmosphere, the telescope, and the instruments on our observed 
contrast measurements. To accomplish this we employ a speckle reconstruction technique espe- 
cially adapted to provide a realistic reconstruction of image amplitudes, taking into account the 
impact of the telescope's Adaptive Optics (AO) system on the Speckle Transfer Function (STF). 
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Employing the G band for a comparison between synthetic and observed contrasts has several 
advantages. Imaging in the band is commonly used in high resolution observations and most ob- 
servatories nowadays have filters with similar passbands, so that comparisons can be readily made. 
Secondly, because a typical G-band filter is 1 nm wide (FWHM), it integrates over many lines, 
masking small uncertainties in individual line positions and averaging out random uncertainties in 
line strengths. The integrated G-band intensity is also insensitive to Doppler shifts and Zeeman 
splitting, allowing a comparison of intensity contrast only. The CH molecular lines that dominate 
the G band form in the photosphere where the approximation of Local Thermodynamic Equilib- 
rium (LTE) applies, which simplifies their numerical modeling. Moreove r, the response of the in- 
tegrat ed G-band intensity to temperature is dominated by the continuum (|Uitenbroek & Tritschlei 
20061), to which LTE definitively applies. 



In their overview of realistic solar convection simulations IStein & Nordlundl (|2000h compare 
histograms of granular contrast observed with the 0.5 m Swedish telescope in the continuum 
with the contrast from their simulations, and find very good agreement after smearing the arti- 
ficial data with a Point Spr ead Function (PSF) representative of the telescope and atmosphere. 
Asensio Ramos et al.l (|2006r) compare the statistics of spatially resolved line profiles of the Fe I 
709.0 nm line observed with the Interferometric BIdimensional Spectrometer (IBIS) at the DST 
with those of profiles computed through a three-dimensional hydrodynamic snapshot, and find very 
good agreement for the intensity distributions in the continuum (after the theoretical profiles were 
adjusted for degradation by the telescope and earth atmosphere) and for the correlation between 
wavelength dependent intensity and intensity in the continuum. However, they also find that the 
spread in velocities in the simul ation is highe r than in the observed profiles, with the latter lacking 



velocities beyond ±1.2 km s . iKellen (|2007h provides an overview of the promises and pitfalls of 



comparing high spatial resolution observations with magneto-hydrodynamic simulations. 

In §[2] we introduce the observational details and in §[3]the speckle reconstruction method. The 
simulations are described in §|4l In §[5] we compare the theoretical and observed G-band contrast. 
Conclusions are given in § [6l 



2. Observations 



The observations were obtained on October 24, 2005 at the Dunn Solar Telescope (DST) 
of the National Solar Observ atory at Sacrarn ento Peak, Sunspot, taking advantage of the high- 
order adaptive optics system (IRimmeld 120041) and a high-speed 2kx2k detector with 12 /i pixels 
manufactured by DALSA for data acquisition. The optical setup provides a field-of-view (FOV) of 
about 61"x61". To select the G-band wavelength range a l.Onm wide interference filter centered 
around 430.5 nm located in front of the detector was used. The diffraction limit (A /D) of the DST 
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Fig. 1 . — Combination of the reference frame {left half) and burst average {right half) of burst 
No. 9. 

at this wavelength is Of'l 17 which together with the detector image scale of Of'OSOpixel^^ leads to 
an approximate twofold spatial oversampling of four resolution elements of X/D per pixel, where 
D — 0.76 m is the aperture of the DST. 



3. Speckle reconstruction 



In order to allow for speckle reconstruction the data was recorded in bursts of 80 images, 
which were acquired within 3.5 sec and streamed off within 30 sec. The exposure time for an 
individual frame was 10 ms. The whole data set comprises 103 bursts with a mean cadence of 
67 sec, covering a time perio d of almost 2 h in total. We employ a speckle masking technique, 
implemented by (IWogei]|2007l) . that incorporates an it erative weighte d least-squares fitting method 
for phase reconstructionand based on a technique by (lMatsonlll99lh . In order to get an improved 
amplitude reconstruction, and thus more reliable photometry, the reconstruction code takes into 
account how the wavefront correction of the AO system affects the STF. 

Each burst frame was divided into 15x15 of l'.'6S square each. The image was then re- 
constructed on each subfield individually, after which the corrected subfields were merged back 
together to give the complete reconstructed frame. A byproduct of the speckle reconstruction code 
is an estimation of the Fried parameter over each of the 15x15 subfields in the FOV. These 
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values characterize the atmospheric seeing over the duration of each burst with higher values cor- 
responding to more favorable conditions. The top panel of Figure [2] shows the mean value of ro for 
each burst averaged over the subfields {solid line and diamonds) and the standard deviation in the 
mean {grey shading). A histogram of ro values per subfield and per burst is shown in the bottom 
panel of Figure [2l Both panels make clear that the atmospheric seeing conditions were good to 
very good during most of the 2 h sequence. 

The top panel of Figure [3] shows the rms contrast variation of a quiet area just below the pore 
(indicated by the white rectangle in Figure [B in all 8240 frames of the 103 bursts. The typical 
uncorrected rms G-band contrast varies between 5 and 7.5% in these individual frames, very rarely 
reaching values over 8%. The selected quiet area is close to the AO lock point, which was set at the 
pore, and is relatively devoid of bright points, and thus is less likely to have strong magnetic fields. 
The speckle reconstruction algorithm was successful in 101 of the bursts. In these successfully 
reconstructed images the rms contrast over the same region of quiet Sun is plotted in the bottom 
panel of Figure [3] (Ja^/zeJ line with open circles). An alternative estimate for the quiet-Sun rms 
contrast, more representative of the contrast in the who le frame and not only close to the lock 



point, is that provided by a blob finding algorithm (see iTritschler & Uitenbroekl 120061) . We use 
this algorithm to locate G-band bright points and eliminate the pixels encompassed by them from 
the rms estimate. The resulting rms values are plotted with the solid curve and filled circles in the 
bottom panel of Figure [3l In general the two estimates for the quiet-Sun rms contrast match each 
other very well. 

Typically, the reconstruction process raises the rms contrasts in the images by about a factor 
of two, to typical values of between 7 and 14%, with the highest contrast over the whole image 
as recovered with the blob algorithm at values close to this upper bound. The contrast of the re- 
constructed images is highly correlated with the contrast of the burst images. During "bad seeing" 
(as judged by the correspondingly low values of ro in Figure [2l top panel) high spatial frequency 
information is lost and cannot be recovered by the reconstruction algorithm. Otherwise, the recon- 
structed contrast would be consistently close to 14%. The upper limit is reached for a range of 
"input image contrasts" of 7-9%. This means the reconstructed contrast saturates, which can only 
be expected if the spatial power spectrum of the granulation does not have a high frequency tail at 
spatial scales below the resolution limit of the telescope. 

With the doubling of the rms contrast, the improvement in image quality is clearly illustrated 
by the sample images in Figure HI The figure shows, from right to left, the central part of a 
frame representing the average over the 80 exposures in one burst (number 15, with the highest 
average contrast in individual exposures and in the reconstructed image of the burst), the reference 
frame for the reconstruction (the image with the highest contrast within the burst), and finally the 
reconstructed image. 
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Fig. 2. — The average Fried parameter ro over the FOV, as determined for each burst by the speckle 
reconstruction code is shown in the top panel, together with its standard deviation (gray shading). 
The bottom panel shows the histogram of tq values determined by the speckle code for all 15 x 15 
subfields in all bursts. 
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Fig. 3. — Rms contrast variation in the individual frames for all 103 bursts {top panel). Vertical 
gray lines indicate the start of each successive burst. Bottom panel shows the quiet-Sun rms con- 
trast of all reconstructed frames, measured from the contrast in a subarea (indicated by the white 
rectangle in Figure [H dashed curve and open circles), and from the blob algorithm (see text, solid 
curve and filled circles). The maximum rms (14.1%) in the reconstructed images derived with the 
latter is indicated by the horizontal dashed line. The vertical dotted lines indicate the two bursts 
for which the iterative reconstruction failed to converge (burst 50 and 72). 
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Fig. 4. — Average (right panel), reference (middle), and reconstructed (left) images of the central 
part of the field-of-view centered on the adaptive optics lock point of burst No. 15. The curve at 
the bottom of each image represents the G-band intensity in a cross section of the image through 
the pore aty = 9.9 arcsec. It illustrates the recovery of the amplitude of high spatial frequencies in 
the reconstructed images. 



4. Simulations 



We compare the speckle reconstructed G-band images with intensities obtained from two dif- 
ferent simulation snapshots that span a realistic range in average field strengths on the solar surface 
outside sunspots. The first one, here referred to as quiet-Sun has an average vertical magnetic field 
of 30 G, which is probabl y at the lower limit of what can be expected in real quiet-Sun areas 
(ITrujillo Bueno et al.ll2004f) . The second snapshot is representative of solar plage with an average 
vertical field of 250 G. It will be referre d to as the plage snapshot . Both snapshots were taken from 
magneto-hydrodynamic simulations by IStein & Nordlundl (|2006b . We empl oy carbon and oxygen 
abundances of log£c = 8.39 and log Co = 8.66, respectively as advocated by lAsplund et al.l (l2005b 
on the basis of C I, C H, and C? line r nodeli ng in three-dimensional hydrodynar nic simulations, an d 
similar modeling by lAsplund et al.l (|2004h of O I and OH lines, and CO lines IScott et al.l (|2006r ). 
The precise value of the oxygen is relevant for the opacity in the G band region because most 
carbon in the cool parts of the solar atmosphere is bound in CO molecules limiting the amount 
available for CH. 

A description of the spectral synthesis is given in ( Uitenbroek & Tritschleijl2006 ). For each of 
the two snapshots the emergent spectrum in the G-band region was calculated at 600 wavelengths 
equidistantly spaced over a range of 3 nm centered at 430.5 nm. In the computation this wavelength 
interval contains 356 atomic lines and 424 electronic molecular transitions of the CH A^A-X^Il 
system with vibrational quantum numbers v — v' = — 0, 1 — 1, and 2 — 2. The calculated spectra 
were then integrated over a typical G-band filter function of 1 nm FWHM with the shape of a 
second order generalized Lorentzian, representative of a dual-cavity interference filter. The rms 
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Fig. 5. — Synthetic G-band filtergrams from quiet-Sun {left panel) and plage (right) snapshots. 
Indicated in white are the rms contrast variations over the whole filtergrams, and indicated in grey 
lettering the rms contrast of pixels outside bright points in the plage snapshot. 



intensity differences of the spatially averaged spectrum from the snapshots with the observed av 



erage solar spectrum are only 7% over the 3 nm spectral range (|Uitenbroek & Tritschlerll2006f) . 
confirming the realism of the synthetic spectra. The main differences between computed and ob- 
served spectra result from slight differences in the central wavelengths of the included spectral 
lines. These intensity differences are inconsequential to the wavelength integrated G-band signal. 

Synthetic filtergrams of the two simulation snapshots as well as the values of their respective 
rms contrast in the G band are shown in Figure [51 The plage snapshot (right panel) includes con- 
centrations of magnetic field in its inter- granular lanes that are strong enough to give rise to G-band 
bright points, while this not the case anywhere in the quiet-Sun snapshot, despite the presence of a 
30 G average vertical field. As a result of the bright points the range in co ntrasts in the plage snap 



shot i s much higher, and the overall rms contrast is raised as well (see also lTritschler & Uitenbroek 



2006L table 2 therein). However, the contrast of the plage snapshot measured in the granulation 
alone, after eliminating pixels in bright points, is much reduced compared to the quiet Sun due to 
the presence of the stronger magnetic field, which inhibits convective flows, and hence heat advec- 
tion. The granular rms contrast is only 16.3% compared to 21.5% in the quiet-Sun snapshot (see 
annotations in grey and white, respectively in Figure [5]), keeping in mind that the former number 
is a slight underestimate because the inter-granular lanes are underrepresented in the rms through 
the elimination of bright-point pixels, which occur preferentially in the dark lanes. 



5. Comparison of observed and theoretical G-band contrasts 




contrast [%] 

Fig. 6. — Histograms of G-band contrast values in the synthetic quiet-Sun filtergram (grey curve) 
and best reconstructed image (No. 15, black curve). 

Comparing the rms contrast values between the synthetic filtergrams and the values for the 
reconstructed images (Figure [3l bottom panel) the large differences are obvious. Where the simu- 
lations have a granular rms contrasts of 21.5% for quiet Sun, the highest contrast in the most quiet 
regions of the observed reconstructed G-band images is 14.1% This difference is further exempli- 
fied in the comparison of histograms of contrast values in the synthetic quiet-Sun filtergram {black 
curve) and the reconstructed frame with the highest contrast (from burst No. 15, gray curve). The 
reconstructed image lacks both relatively dark pixels with contrasts less than —30%, and relatively 
bright pixels, with excess contrast of +30% compared to the synthetic filtergrams. In addition, the 
histogram of the synthetic filtergram is more clearly a bimodal distribution of two Gaussians with 
a main peak at —20% and a secondary peak at -|-25%. The histogram of the observed filtergram 
can also be closely represented by a sum of two Gaussian distributions, but with maxima that are 
much closer zero, giving the impression of a unimodal rather than a bimodal distribution. 

The darkest (contrast less than —0.30, white contours) and brightest features (contrast less 
than -1-0.30 black dashed contours) in the synthetic quiet-Sun snapshot are outlined in Figure |7l 
Clearly, the dimensions of the darkest and brightest areas in particular are well above the theoretical 
resolution limit of the DST (which is O'.'lll at the central wavelength of the G band), and should 
be resolvable in the speckle reconstructed images when seeing is good. This suggests that the 
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Fig. 7. — Contours outline the pixels with contrasts of less than —30% (white contours), and more 
than +30% {black contours). 



properties of the simulation we designate as quiet-Sun are inherently different from the real quiet 
Sun. In turn this implies that temperature fluctuations in the simulations are overestimated because 
granular intensities in the G-band filter signal are do minated by the continuum and are directly 



related to temperature at heights between 25 - 50 km. dUitenbroek & Tritschlerll2006h . 



To verify that the quiet-Sun simulation snapshot does not represent an exceptional case of rms 
high contrast, we investigated the statistics of rms contrast of subfields of comparable size in the 
observations. A histogram of rms contrast values in 8!'5 x 8"5 subfields of the 101 successfully 
speckle-reconstructed images is shown in Figure [8l Only 0.89% of the subfields have a rms con- 
trast of 14.1% (the highest value of rms in the quiet area of the reconstructed frames), and only 
0.01% has an rms of 20% or more. It is therefore highly unlikely that the simulation snapshot 
we employed happens to correspond to the high contrast tail of the distribution of subfield rms 
values. At first sight it seems questionable that less than 1% of subfields have rms contrast more 
than 14.1%, since there are several frames with quiet-Sun rms contrast that high (see Figure [3l 
bottom panel). The explanation for this seeming discrepancy is the 5-minute oscillation, which has 
patches of intensity variation of typical size 10"-15". Such large patches contribute much more to 
the overall rms in the frames than to the rms contrast in subfields smaller than 10". 
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Fig. 8. — Histogram of rms contrast variation over 8.5x8.5 subfields in the quiet-Sun are of the 
reconstructed images. 



One type of image degradation that reconstruction techniques do not have a handle on is 
the effect of spatially unresolved scattering, essentially the image smearing that occurs because 
of the broad outer wings of the point spread function (PSF) of the telescope and atmosphere. Is 
this perhaps the main reason for the substantially lower contrast of the reconstructed images with 
respect to the synthetic filtergrams? Below we present a simple argument that this is very unlikely 
the case. 

Let lij the intensity in a synthetic solar image and (/) = j^Ylij^ij its average over the A/^ pixels 
of the image. The rms contrast C of the intensity is defined as: 



In the simplest approximation scattering takes a fraction a of the spatially resolved intensity and 
spreads it equally over all pixels. The resulting degraded image /' can therefore be written as 
/' = (1 — «)/ + «(/), with (/') = (/). The rms contrast of the degraded image is 



5.1. Spatially unresolved scattered light 




(1) 




(2) 



For the calculated contrast to be degraded from the nominal 21.5% to the observed contrast in the 
best reconstructed images of 14.1% requires a scattering coefficient a of at least 0.34. This is 
unrealistically high for the DST by at least a factor of three, as is judged from observations off the 
solar limb, and a transit of the planet Mercury (Reardon, private communication). 



6. Conclusions 

We have compared the rms intensity contrast in the G-band in a series of speckle reconstructed 
images obtained at the D = 0.76 cm DST with that in two synthetic G-band filtergrams computed 
from magneto-hydrodynamic snapshots representative of quiet-Sun and plage with average vertical 
magnetic fields of 30 G and 250 G, respectively. We find that the rms G-band contrast in the best 
speckle-reconstructed images is just over 14%, as judged from a quiet subarea, relatively devoid 
of bright points, and from the overall contrast in the image after removing bright points (Figure [3l 
lower panel). This observed rms contrast is contrast is considerably lower, by more than 50%, than 
the rms contrast in the synthetic quiet-Sun filtergram which was found to be 21.5%. In the plage 
snapshot the rms contrast of the granulation, excluding pixels with bright points, is 16.3%. 

The employed speckle reconstruction algorithm accounts for the effect of the telescope's AO 
system, and provides reliable photometry in the reconstructions. During moments of consistently 
good seeing over the duration of an 80 frame speckle burst (corresponding to typical values of the 
Fried parameter ro of 12 cm or larger, compare Figures |2] and [3l top panel), the reconstructed rms 
contrast seems to saturate to 14.1%, suggesting this approaches the contrast on the real Sun. To 
verify the conclusion that the rms G-band in the observed quiet Sun is indeed substantially dif- 
ferent from the that in the synthetic snapshot we convolved the latter first with the telescope Point 
Spread Function (PSF) appropriate for the D = 0.76 m DST, and then with a PSF representative for 
telescope plus ro = 15 and ro = 10 cm seeing. The resulting rms contrasts were 19.3%, 16.8%, and 
14.9%, respectively. These values are again more than 50% higher than the rms contrast in the best 
quiet-Sun subarea in the uncorrected single frames (Figure |3l top panel), which do not reach more 
than 9%, even though they are estimated close to the AO lock point. Remarkably, the rms contrast 
of the quiet-Sun synthetic filtergram, degraded to the values representative of far from superior 
seeing (i.e., ro = 10 cm), is higher than that in the best reconstructed images, which should have 
all image degradation, apart from spatially unresolved scattering, taken out. 

Areas in the synthetic filtergram with negative contrast in excess of that in the observed quiet 
Sun are found in the inter-granular lanes, those with positive excess in relatively large (several 
tenths of arcsec) patches in the granules (see Figure |7]). The spatial scales of these areas are 
large enough, for positive excess in particular, so that they should be resolved in the best speckle 
reconstructed images. The synthetic filtergram, therefore, does not show significant power that 
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contributes to the rms at spatial scales smaller than the resolution limit of the DST, which is 0''1 17 
in the G band. The saturation of the reconstructed rms contrast to just above 14% during the mo- 
ments of best seeing seems to confirm that no significant power is present in the solar granulation 
at spatial frequencies close to the resolution limit of the telescope. The simple consideration of 
degradation of contrast by spatially unresolved scattered light presented in Section [STI shows that 
it is very unlikely the cause of the discrepancy between computed and observed G-band contrasts. 
We remark that, given the typical spatial scale of intensity fluctuations caused by the 5 min oscilla- 
tions, the latter contribute more to the rms of the large FOV spanned by our images than to the rms 
contrast in small fields like our simulation snapshots, which are smaller than those p-mode spatial 
scales. Therefore, the discrepancy in G-band contrast caused by convective motions alone is even 
larger than it seems from our estimates. We conclude that there remains a clear and unresolved 
difference between the simulated and observed rms G-band contrasts: the quiet-Sun simulation 
snapshot overestimates the rms G-band contrast by as much as 50%. This raises the fundamental 
question whether the simulations still lack fundamental physics or sufficient numerical resolution, 
or whether we do not fully understand the degradation of contrast by the earth atmosphere and 
telescope during observations. 

One possible source of the contrast discrepancy is the absence of cooling by CO molecules 
in the simulation, where radiative transfer is approximated by only four op acity bins, and four 



direc tions. CO is the dominant coolant above the photospheres of cool stars lAyres & Testerman 



(e.g., 119811) . and its presence may well reduce temperature differences in that layer, leading to a 
reduction in rms contrast variation in the G band since the latter is directly related to temperature 
variations at about 25-50 km above the photosphere. Ho wever, the tempe rature reduction by CO 



may only be important in layers that are higher than that dScott et al.ll2006h . 



The clearest hint, however, towards a solution to the G-band contrast discrepancy comes from 
the much reduced granular contrast in the plage simulation snapshot. Clearly, the presence of 
high concentrations of magnetic field inhibits convective motions and convective heat transport, 
and reduces tempera t ure di fferences and thereby, G-band contrast. The findings described by 



Asensio Ramos et al.l (120061) that s i mulate d profiles of the Fe I 709.0 nm line in a purely hydro- 
dynamic simulation I Asplund et al. ( 2000a ) show high velocties in excess of ±1.2 km s^^ that are 



lacking in high resolution observed quiet-Sun observations confirms that such simulations without 
magnetic field seem to overestimate convective flow velocities. We speculate that the real quiet 
Sun contains much more magnetic field than the average 30 G that we assumed in our "quiet-Sun" 
snapshot, so that the comparison of its rms G-band contrast to that in the observed quiet-Sun is not 
appropriate, much like comparing apples to oranges. This magnetic field has to be present mainly 
in the form of weak fields, since it does not show up in the form of bright point in the G band. We 
further speculate that, if indeed the real Sun contains more magnetic field and if convective flow 
speeds are less than what is derived from field-free simulations, the element abundances that are 
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derived from such simulations are underestimated. For slower convective flows abundances would 
have to be increased to match the observed widths of lines. This would alleviate, at least in part, the 



conflict that exis ts between the low abundances of oxygen and carbon (|Asplund et al.ll2005L l2004l : 



Scott et al]|2006 ). derived from three-dimensional simulations compared to their meteo ritic values. 



and th e conflict these low values pose with current models of the solar interior (eg., lAyres et al. 



20061). 



For more detailed comparison between observed and simulated G-band contrast more simu- 
lations are required with intermediate field strengths between the 30 G and 250 G average values 
we employed here. In addition, observations from space, where the degradation of image quality 
by the earth atmosphere is eliminated, and only the easier to characterize contrast reduction caused 
by the telescope and instrument has to be considered, should conclusively establish the value of 
the G-band contrast of the real Sun, hopefully confirming our conclusions here. This has to be 
accompanied by deep measurements of the magnetic (vector) field, so that we can finally answer 
the question what really constitutes quiet Sun. 



We are grateful to Bob Stein for providing the simulation snapshots, and to Friedrich Woger 
for allowing us to use his speckle reconstruction code. 
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